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We study a singlet-triplet Josephson junction between a conventional s-wave superconductor and 
an unconventional px-wave superconductor. The Andreev spectrum of the junction yields a spon¬ 
taneous magnetization in equilibrium. This allows manipulating the occupation of the Andreev 
levels using an ac Zeeman field. The induced Rabi oscillations manifest themselves as a resonance 
in the current-phase relation. For a circularly polarized magnetic field, we find a spin selection rule, 
yielding Rabi oscillations only in a certain interval of the superconducting phase difference. 


The current-phase relation of a Josephson junction 
contains information about the Andreev levels and their 
occupations. Junctions formed between unconventional 
superconductors have exotic bound states leading to un¬ 
usual current-phase relations. Amongst unconventional 
Josephson junctions, those realized between singlet and 
triplet superconductors are of special interest, because of 
their incompatible spin pairing symmetries. Their equi¬ 
librium properties have been studied for various types of 
heterogeneous junctions [IH3. 

Let us consider such a Josephson junction between 
a conventional spin-singlet, s-wave superconductor and 
an unconventional spin-triplet, p 3 ,-wave superconductor. 
This junction displays exotic spin properties. Namely, 
it hosts two spin-polarized Andreev bound states, which 
have the same spin Q. In equilibrium, this results in 
a spontaneous magnetization that is 27r-periodic in the 
superconducting phase difference. On the other hand, 
a TT-periodic equilibrium supercurrent, which does not 
probe the exotic spin properties, is found The spin 

properties of the Andreev levels open the possibility for 
spin manipulation, using a time-dependent Zeeman field. 
A similar idea, the manipulation of the Andreev levels 
in spin active Josephson junctions between conventional 
superconductors, has already been reported Uil- 

In this article, we show, that an ac Zeeman field leads 
to coherent Rabi oscillations between different spin states 
of the singlet-triplet junction. These Rabi oscillations 
manifest themselves as resonances in the current-phase 
relation. For a circularly polarized magnetic field, we 
find a spin selection rule, yielding Rabi oscillations only 
in a certain interval of the superconducting phase dif¬ 
ference. The applied Zeeman field also induces non¬ 
coherent transitions between the bound states and the 
continuum states. In principle, these transitions, which 
we treat within a macroscopic master equation approach, 
could give rise to a decay mechanism for the Rabi oscil¬ 
lations. However, we find, that due to spin and energy 
constraints, these processes do not coexist with the Rabi 
oscillations. 

A possible experimental realization of our proposal 
could be based on the quasi-onedimensional (TMTSF) 2 X 
Bechgaard salts [l^ in which a pa;-pairing was discussed. 



FIG. 1. Setup of an effective singlet/triplet junction using 
a semiconducting ferromagnetic nanowire (F NW) contacted 
with conventional singlet superconductors (S). The gate (G) 
allows for realizing a barrier with tunable transparency. By 
adjusting the length of the nanowire between the gate and 
the right superconductor, one can realize an effective triplet 
superconducting reservoir. 


as suggested in Ref. Q- Alternatively, we propose to 
realize a junction between conventional superconductors 
separated by a ferromagnetic semiconducting nanowire, 
as illustrated in Fig.[TJ The gate would allow for the real¬ 
ization of a barrier with tunable transparency. Further¬ 
more, an effective pa;-wave superconductivity is realized 
when the length of the nanowire between the gate and 
one of the leads matches the coherence length for the 
superconducting correlations induced in the nanowire. 

Let us introduce the model. The Hamiltonian describ¬ 
ing a Josephson junction between an s-wave supercon¬ 
ductor and a one-dimensional, time-reversal symmetric 
Px-wave superconductor reads 

= ( 1 ) 

where = (Rj,, l|, Lj,, R|)^, and Rl and Lj. are creation 
operators for right-moving and left-moving electrons with 
spin a respectively. The Bogoliubov-de Gennes 

Hamiltonian Pi is given as 

Pi = VfPVzTz +U{x)riz;Tz - As(x)tx -b , 

( 2 ) 

where Tx^y^z and r]x,y,z denote Pauli matrices in particle- 
hole and R/L spaces, respectively. The first term in 
Eq. (I2|), with Fermi velocity wp and momentum opera¬ 
tor p, is the kinetic energy. The second term describes 
normal backscattering induced by a scalar potential U (x) 
that vanishes outside the central region of the junction. 
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0 < X < L, where L is the junction length. It gives rise to 
an electronic transmission probability T, when the junc¬ 
tion is in the normal state. The third term describes 
s-wave pairing with gap As(x) = AsO(—x), where 9 is 
the Heavside step function, on the left side of the junc¬ 
tion. The last term describes Px-wave pairing with gap 
Ap(x) = Ap9{x — L) between electrons having opposite 
spins along z-direction on the right side of the junction. 
Without loss of generality, Ag and Ap are assumed to 
be real. For simplicity, we will restrict our analysis to 
the case Ag = Ap = A [l3|. The superconducting phase 
difference across the junction is denoted (j). Note that we 
use units where h = 1. 

In the short-junction limit, L <C r'p/A, we find that 
the Bogoliubov-de Gennes Hamiltonian (ED admits for 
two bound states with energies 0 ] 
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FIG. 2. Energy-phase relation of the two bound states for the 
transmission T = 0.8. The continuum of states is indicated 
in gray. All states shown have spin f. 


, (3) 

and wavefunctions '4’±(x) that are given in the Supple¬ 
mental Material [l^. Note that the choice of the spinor 
df implies that we are considering states with spin down 
only. Furthermore, Eq. (ED admits for a four-fold degen¬ 
erate continuum of (outgoing) propagating states with 
energies E {\E\ > A) and wavefunctions ijjEtiix), where 
p. is a degeneracy index. Using a Bogoliubov transforma¬ 
tion. 


= + ^'pEt,{x)'yEti , (4) 

I'—zt E,fi 


where and 7 ^;^ are annihilation operators for quasi¬ 
particles in the bound state with energy E^ and for quasi¬ 
particles in the continuum with energy E and degeneracy 
index /x, respectively, we may diagonalize the Hamilto¬ 
nian (ED to obtain 




( 5 ) 


u=± 


E ,(-i 


A typical spectrum is shown in Fig. El Note that at van¬ 
ishing coupling, T —>• 0, the spectrum of the s-wave lead is 
gapped, while the Px-wave lead, which realizes two copies 
of the Kitaev model 1^ in opposite spin sectors, admits 
for a zero-energy edge state. A finite coupling moves this 
state to finite energy and yields the bound state e = —, 
while a second bound state {e = -b) detaches from the 
continuum. In contrast to conventional junctions, both 
bound states carry the same spin (ct = 4 .). 

The bound state occupations, = {'ylju), determine 
the magnetization carried by the junction. 


M = - 


fJ-B 

2 



where /xb is the Bohr magneton. In equilibrium, = 
f{E^), where / is the Fermi function. As a result, the 


junction carries a spontaneous magnetization, which is 
27r-periodic in the phase difference [7|. 

The Josephson current is given as 



The equilibrium supercurrent is spin-insensitive and, 
therefore, it is xr-periodic. 

Thus, to probe the peculiar spin properties of the junc¬ 
tion, we have to consider out-of-equilibrium effects. In 
particular, in order to manipulate the bound state occu¬ 
pations, we apply an ac magnetic field, which is described 
by the Zeeman Hamiltonian 

Hz = hb Y 

. 

We assume the magnetic field to be uniform and circu¬ 
larly polarized in the plane perpendicular to the spin 
quantization axis, B = H(cosHt, sin 17^,0), where H is 
the driving frequency. Such a field leads to spin-flip pro¬ 
cesses. In particular, the spin of the system changes by 
ASz = sgn (H) when a photon is absorbed, whereas it 
changes by ASz = —sgn(H) when a photon is emitted. 
In the following, we concentrate on the case H < 0. 

In order to identify the processes induced by the field, 
we express Eq. (ED in terms of the quasiparticle operators 
using the Bogoliubov transformation (ED- We find 


J dx B ■ [rIcTss'Rs' + LIcTss'Lsi) 


Hz = E+,_7+7- + E 




+ ^ E I + h.c. , (9) 


E,E'\fi,fi' 


where 14 ,a' = J dx tjj\r]z;{-iTy)tjj\' for A S {-b, -, Efi}. 

According to Eq. m, the field couples two quasipar¬ 
ticle states. Three different types of processes are possi¬ 
ble: transitions involving only bound states [first term of 
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FIG. 3. a) The four possible states of the junction and 
their magnetization M. Full (open) dots represent occupied 
(empty) states, b) Transitions induced by a circularly polar¬ 
ized magnetic field with < 0. The shaded region is the 
continuum of states. The thick black arrows denote the spin 
of the states. Absorption (emission) of a photon changes the 
spin by ASz = —1 (ASz = -fl). The Andreev bound states 
carry a spin down (up) in the phase interval 0 < (j) < tt 
{■n < (f) < 27r). 


Eq. dSD], transitions involving a bound state and a con¬ 
tinuum state [second term of Eq. ([5])], and transitions 
involving only continuum states [third term of Eq. ([^]. 

For the discussion of the spin properties, we note that 
the destruction of a quasiparticle with spin down at nega¬ 
tive energies corresponds to the creation of a quasiparticle 
with spin up at positive energies. So far, we used both 
positive and negative energies for spin j,. In the following, 
we will work with both spin directions, but only positive 
quasiparticle energies. Furthermore, we will assume that 
the temperature is low (on the scale of A), such that the 
continuum states at E > A are empty. 

The transitions involving only bound states correspond 
to Rabi oscillations, i.e., coherent oscillations between 
the state |0), where both bound states are empty, and 
the state |2), where both bound states are occupied. 
Rabi oscillations occur when the oscillation frequency |n| 
matches the Rabi frequency, flu = -\-E-{(j))\. Us¬ 

ing the energy dispersions given by Eq. the inter¬ 
esting frequency regime is, thus, given by A < |r2| < 
V2A 0. Then, sweeping the phase at fixed frequency, 
the resonance condition is met for 4 different values of 
the phase: ())o, tt — (/)o, tt -|- (/iq, 27r — ())o. However, the spin 
selection rule imposes a further constraint. Namely, at 
U < 0, Rabi oscillations are possible only if the bound 
states carry a spin down, which is the case in the interval 
Q < (j) < TT. Thus, the circularly polarized magnetic field 
leads to Rabi oscillations only at two values of the phase: 
4>q and TT — (j)Q. 

Transitions involving a bound state and a continuum 
state change the parity of the bound state occupation. 


namely they connect the even-parity subspace {|0), |2)} 
to the odd-parity subspace {|—), |+)}, where \v) denotes 
the state in which only the bound state with energy \Ei,\ 
is occupied. (For a sketch of the four states of the junc¬ 
tion, see Fig. [3^).) This would represent a decay mech¬ 
anism for the Rabi oscillations. We may distinguish two 
different processes, that are sketched in Fig. [3 ]d). 

In an ionization process, a particle from a bound state 
is promoted to a continuum state. Energy conservation 
imposes the condition |fl| > A — \Ev\ for such a process. 
In the frequency range of interest for Rabi oscillations, 
this condition is always met. However, the spin selection 
rule imposes that, at H < 0, the bound state carries a 
spin up, which is the case in the interval tt < (p < 2 tt 
only. Thus, Rabi oscillations and ionization processes 
occur in different phase intervals. 

In a refill process, a Cooper pair is broken such that 
one particle occupies a bound state, whereas the second 
particle is promoted to a continuum state. Here energy 
conservation imposes the condition |H| > A -|- In 

the frequency range of interest for Rabi oscillations, this 
condition is never met for the state with energy iF+j. 
By contrast, for the state with energy |F_|, one obtains 
a critical phase pc such that the condition is met in the 
phase intervals [—pc^Pc] and [tt — pc, tt + pd- Compar¬ 
ing pc with the phase po giving the resonance condition 
for Rabi oscillations, we find pc < po- Thus, Rabi oscil¬ 
lations and refill processes also occur in different phase 
intervals. Note that, here, the spin selection rule imposes 
that, at H < 0, the bound state carries a spin down, which 
is the case in the interval 0 < p < tt. 

We conclude that transitions between a bound state 
and the continuum due to a circularly polarized magnetic 
field do not provide a decay mechanism for the Rabi os¬ 
cillations due to energy and spin constraints. However, 
such a decay may be due to other parity non-conserving 
processes related to, e.g., quantum phase fluctuations due 
to the resistive environment of the junction [13,0. 

Finally, transitions involving only continuum states 
have a threshold |H| > 2A. Thus, they do not play a 
role in the frequency range of interest for Rabi oscilla¬ 
tions. 

The modifications of the bound state occupations, 
which are induced by the different processes discussed 
above, may lead to strong deviations of the Josephson 
current o from its equilibrium value. To compute the 
Josephson current in the steady state, we introduce the 
matrix elements pap = {a\p\j3) of the reduced density 
matrix p, where \a),\j5) G {|0), |-|-), |—), |2)}, so that 
Eq. Q reads 


/ = {I++I-){poo-p22) + {I+-I-){.P—-p++) , (10) 
where I± = —e{d\E±\/dp). 

Taking into account the above considerations for the 
frequency range where Rabi oscillations can take place. 
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we find the following behavior in different phase inter¬ 
vals. For phases (j) G [tt, 27r], only ionization processes 
are possible. Thus, the bound states are always empty, 
i.e., Poo = 1 and p _= p++ = P 22 = 0. As a con¬ 

sequence, the current remains equal to its equilibrium 
value, -\-1_. On the other hand, in the phase 

intervals [0,0c] and [tt — 0 c, tJ"], the ac field yields a re¬ 
fill process for the state v = —. Assuming that the 


rates for parity non-conserving processes due to the en¬ 
vironment are much smaller than the field-induced rates, 

this state then will always be filled, i.e., p _= 1 and 

Poo = P++ = P 22 = 0. Thus, the current is / = 1+ — I-. 

To evaluate the current in the phase interval [0c, tt — 
0c], which includes the phases 0o and tt — 0o where Rabi 
oscillations take place, we use the stationary solution of 
the master equation. 
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Here, wi = 2V+ _pbH with = r^|sin0|(l -|- 

|sin0|)/(l -I- r |sin0|)^, 5uj = O-|-sgn(sin0)OR,, po 2 = 
e*^Vo 2 , and p 20 = e“*^*p 20 - Furthermore, f]/^ are the 
ionization (I) and refill (R) rates of the state v, respec¬ 
tively, and Fs = +1"^ )• The other 10 elements 

of the 4x4 density matrix that are not shown remain 
zero along the time-evolution. 

Sol ving Eq. (ITT1) to obtain the steady-state occupa¬ 


tions 


14| . we find that the current. 


I = r 


r2 -h (2,5a;)- 




( 12 ) 


is the sum of a background term and a reso¬ 
nant term. The background term is given as = 
E.=± - r]})/r., where F± = F^ + Fg. The cur¬ 

rent at resonance is given as 


cO ^ 

r 2 


rs/° 


1 ^ 11 ^ 

r+r_ 


(/+ -i-)Y. - r' 


u=± 


(13) 


whereas the width of the resonance is determined by 


As pointed out above, all field-induced decay rates are 
zero in the phase interval [0c, tt — 0c]. Therefore, we 
introduce phenomenological rates 7 to describe the parity 
non-conserving processes due to the environment. At low 
temperature, the refill processes are negligible [l^, and 
we are left with two rates 7 J,. In that case, I°° reduces to 
the equilibrium current, I°° = Assuming 7 ^ <C |wi|, 
the current at resonance is obtained as 

( 15 ) 

7+ +7i 


whereas the width of the resonance is 


F 


7^+7! 



|wi| . 


(16) 


As the state -I- is closer to the continuum, we expect 7 ]j_ > 
7 L. Depending on their relative magnitude, the current 
may be completely suppressed at resonance, when 7 L = 
7 ]|_, or change its sign as compared to the equilibrium 
current, reaching a magnitude ~ /+-/_, when 7 ]_ <C 

7l- 

Fig. SI shows the non-equilibrium current-phase rela¬ 
tion for a circularly polarized Zeeman field in the en¬ 
tire phase range. The 27r-periodicity is due to the spin- 
sensitive manipulation of the bound state occupations. 
If the sign of D was reversed, the current-phase relation 
would be phase-shifted by TT, i.e., /(D, 0) = /(—fl, 0-|-7r). 

By contrast, the spin sensitivity is lost, if we use a lin¬ 
early polarized field, B = 2B{cosV,t,0,0). Such a field 
can be viewed as the superposition of two circularly po¬ 
larized fields with opposite helicities. Thus, there is no 
spin selection rule anymore, and Rabi oscillations may 
now occur at 4 phases (0o, tt — 0o, tt + 0o, 27r — 0o). Fur¬ 
thermore, the field-induced ionization rates are non-zero 
for all superconducting phase differences while F 5 exists 
in the phase intervals [— 0 c, 0 c] and [tt — 0 c, tt + 0 c]. 

Extending the master equation (HH) to the linear case 
and introducing a rotating-wave approximation to de¬ 
scribe the vicinity of the Rabi resonances, we find that 
the steady-state current is given by Eqs. (HHi-dll with 
rf = Fj((D) -h r^(-D) , where X = I, R. 

The transition rates for the ionization and refill pro¬ 
cesses involving the bound state v can be calculated from 
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FIG. 4. Current-phase relation for a junction with transmis¬ 
sion T = 0.8 and a circularly polarized magnetic field, with 
amplitude = 10~^A and frequency Q = —1.3A. The 
phenomenological ionization rates are chosen as ylj. -I- yL = 
10“®A with yL = 0.1y]j_ (solid red line) and yL = y]j_ (dot¬ 
ted red line). The equilibrium current is given for comparison 
(dashed black line). 



— <t> = 3TTl2 

— 4> = b7: f 4. 
- <f> = 9TT/8 


1 I 2 



ini/A 



FIG. 5. Ionization (I) and refill (R) rates for T = 0.8 as a 
function of the driving frequency fl < 0 for several phase dif¬ 
ferences 4>- a) Ionization rate of the state v = b) Ionization 
rate of the state v = —. c) Refill rate of the state v = d) 
Refill rate of the state u = —. 


FIG. 6. Current-phase relation for a junction with transmis¬ 
sion T = 0.8 and a linearly polarized magnetic field, with 
^b.B — 10~^A and = —1.3A (solid red line). The equilib¬ 
rium current is given for comparison (dashed black line). 


is shown in Fig. [6l As the manipulation of the bound 
state occcupations is not spin sensitive in that case, the 
current is 7r-periodic as in equilibrium. Generically, the 
out-of-equilibrium current-phase relation is 27r-periodic 
as soon as the ac field carries a finite angular momentum, 
leading to spin-dependent refill and ionization rates. 

In conclusion, we have shown that the occupations of 
the Andreev levels in a Josephson junction between an s- 
wave and a Px-wave superconductor can be manipulated 
using an ac Zeeman field. The induced Rabi oscillations 
manifest themselves as resonances in the current-phase 
relation. For a given circular polarization, their presence 
or absence depends on the spin state of the junction, thus 
providing a spin detection scheme. 

LE thanks Roman-Pascal Riwar for helpful discus¬ 
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the Universite Grenoble-Alpes, by ANR through grants 
ANR-11-JS04-003-01 and ANR-12-BS04-0016-03, and by 
an EU-FP7 Marie Curie IRG. 


Eq. ([9]) using Fermi’s Golden Rule, 

Fyi^(fl) = /“ dE p{E) |K.y£;sg„(sin0) 

X d[f2-|-(|E,,| =F^^)sgn(sin(/))] . (17) 


Here p{E) = {2ttvp)~^ E j ^E"^ — is the density of 
states in the leads. Fig. 0 shows the frequency depen¬ 
dence of the ionization and refill rates. Note that the 
rates f]/^, whose typical amplitude is ~ van¬ 

ish below the threshold frequency H = A =F |A'!^|) as dis¬ 
cussed above. Furthermore, they are suppressed at large 
frequencies |H| ^ A, while they display a maximum in 
the vicinity of the threshold frequency. Simplified ex¬ 
pression for the rates in the ballistic (T = 1 ) and opaque 
(T = 0) limits are provided in 141. 

The current-phase relation for a linearly polarized field 
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Supplemental Material: Magnetic resonance in a singlet/triplet Josephson junction 

The Supplemental Material provides technical details on some results that were used in the main text. In Sec. HI 
we derive the eigenstates of the Bogoliubov-de Gennes Hamiltonian, Eq. (2) in the main text. In Sec. HIl we use these 
eigenstates to compute the ionization and refill rates due to an ac Zeeman field, Eq. (17) in the main text. In Sec. IIIIl 
we provide the stationary solution of the master equation, Eq. (II) in the main text. 


I. WAVE FUNCTIONS 

In order to obtain the eigenstates of the Hamiltonian (2), we determine the general form of the wave functions 
in the leads in Sec. ITU The wave functions are given in the basis ?7 0 r, where rj denotes the R/L space and r the 
particle-hole space. Then we use the boundary condition at the junction to establish the wave functions for the bound 
states in Sec. m and for the continuum states in Sec. Ilcl We provide simple expressions both in the cases of a 
transparent and an opaque junction. As the wave functions are 27r-periodic, we restrict our considerations to the 
interval (p S [0, 27r[. 


A. Wave functions in the leads 

In the left (s-wave) lead, a; < 0, the Hamiltonian (2) reduces to 

Us = VfPVzTz - (SI) 

It has a block-diagonal structure in the R/L space. In each block, characterized by = ±1, we thus need to solve 
an auxiliary 2x2 eigenvalue problem given by 


{±vfPTz - At^) Q =e(^^ . 


(S2) 


Using the solutions for this problem, we find that the most general form of the wave functions associated with the 
Hamiltonian (fSTI) at energies above the gap, \E\ > A, is the superposition of four independent spinors, 
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(S3) 


Here, a = (E — sig^l{E)^/E^ — A^)/A and k = sigii(E)VE^ — A^/uf- The prefactor in Eq. (IS3I1 ensures that each 
4-spinor is normalized to unity. Furthermore, 
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(S4) 


are nothing but the BCS coherence factors. Thus, the spinors with the coefficients A“, A™*, and A™ in 

Eq. (IS3I) describe right-moving electron-like, left-moving hole-like, left-moving electron-like, and right moving hole- 
like quasiparticles, respectively. 

Below the gap, \E\ < A, there are only two evanescent solutions, such that the most general form of the wave 
functions associated with the Hamiltonian (|S1I) reads 
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(S5) 


Here a = {E — ly/A'^ — E'^)/A, which may be written as a = e where x G R is the phase shift acquired in an 
Andreev reflection process, and k = y/A"^ — E'^ Jvf gives the decay length in the lead. 





























(S6) 


In the right (p-wave) lead, x > L, the Hamiltonian (2) reduces to 

T-Lp = VFprjzTz - Ar]zTxe' 


We notice that 


Tip = U'^'HsU, where U = exp iTz + ^ 2 ) 

This allows us to write the general form of the wave functions both in the continuum. 
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and below the gap, 


ip{x) = 



( ^ 


f 



De 

ae-^^ 
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+ Df 


0 



1 0 J 



1 J 



(S7) 


(S8) 


(S9) 


To determine the coefficients in the wave functions introduced above, we need to match them at the junction. For 
this, we derive the transfer matrix M associated with the scalar potential U{x) = Uq9[x{L — cc)] in the normal part 
of the junction. When Uq is large, the wave functions with energy E in the normal part of the junction, Q < x < L, 
are readily obtained as 





(A 


= E< 

^\e-^^+E: 

Ae^^+E< 

* e-""+U> 



\oJ 

\oJ 

\lj 

KiJ 


(SIO) 


where A = Uq/vf- Using the continuity conditions for the wave functions at a; = 0 and x = L, we can get rid of the 
coefficients , E ^, E ^, E ^, and establish the relation 


^{L) = M^( 0 ), 


(Sll) 


where M = cosh(AL) + sinh(AL)r 7 y. The coefficients in the transfer matrix can be related to the junction transparency, 

T =1! cosh^(AL). (For definiteness, we will assume A > 0 below.) At T = 0, the two superconductors are decoupled. 

In that case, the boundary condition Eq. (ISllI) reduces to 

(l + r;,)^(0) = 0, (SI2) 

(I - 77y)'!/'(T) = 0. (SI3) 

Below we use the matching condition (ISllI) to obtain the bound state and continuum wave functions. Furthermore, we 
consider the short-junction limit, L —> 0, while keeping the product UqL that determines the transparency constant. 


B. Bound state wave functions 

In the transparent case, T = 1, the matching equation provides two solutions in the R and L sectors, respectively. 
The solution in the R sector has energy 


Ur 


A sin (sin 4 >); 


(S14) 


its wave function is obtained with = Be = 0 and Dg = ie*'^'^^sgn (sin(()) Using the normalization condition for 
the wave function, we can fix B^ = \/ A| cos((()/2)|/(2wf) = \/ kr/2. The solution in the L sector has energy 
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its wave function is obtained with = Q and Dh = —where Be = \/A| s\n{(j)/2)\/{2 vy) = -\/kl72- 

The two states cross at </> = •n/2 and <j) = 37r/2. The connection to the energy E+ {E-) given in the main text is 
made by taking for each intervall the state with the higher (lower) absolute value of the energy, i.e., 


£;+(</)) = sgn(sin(())max{|£’R((())| , |Bl((/))|}, (S16) 

E-{(j)) = sgn(sin(;ii)min{|£'R((;ii)| , |Bl((())|}. (S17) 


At finite backscattering, these solutions hybridize and an avoided crossing appears near the phases 4> = ^ and cj) = ^. 
In the opaque case, T = 0, the higher-energy state merges with the continuum while the matching equation provides 
a unique bound state solution with energy E_ = 0 that resides on the right side of the junction only. The coefficients 
are given as Bh = Be = 0 and Bh = e“®‘^Be with Be = /S./{2vy) = \/k_/2. 

At arbitrary transmission, we find two eigenstates with energies given by Eq. (3) in the main text. Using the 
matching condition, Eq. (ISllI) . and the normalization condition, f (ix|'I'(x)| = 1, we obtain the coefficients for the 
bound state with energy E^: 


b: = VTcos (x. + c'^ 

B(( = -^v/T(l-T)cos^C^ 
B3" = sinx.C", 


B 


B _ 

h — 


(l-r)cos^- 


cos 




(518) 

(519) 

(520) 

(521) 


where 




2 cos 2xu{T cos^ I — sim Xu) 


(S22) 


Note, that the expressions previously given for the special cases T = 0 and T = 1 differ by an irrelevant global phase 
factor. 


C. Continuum wave functions 


Eor a hxed energy in the continuum, \E\ > A, the relation between the four incoming and four outgoing wave 
functions encoded in Eqs. (IS^ and (ISSl) can be expressed through a scattering matrix S[E) such that 




i 

\C^^) 


= S-\E) 


^out 

^out 


(S23) 


The scattering matrix is unitary, i.e., S ^ . At energies \E\ ^ A, the scattering matrix simplifies to S' = 

-iVl - Ttz + VTrix, in agreement with the transfer matrix introduced in Eq. (ISllI) . 

Eor a transparent junction, T = 1, the scattering matrix is block diagonal as the R and L sectors decouple. It reads 
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S = 


0 

1 -a" 


a{l+e~'*) 
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l-g" 
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(S24) 


/ 


1 —1 — OL^e^^ 

Eor the opaque junction, T = 0, scattering states are confined within each lead, and the scattering matrix reads 

. / -z 0 0 0 \ 
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(S25) 
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where we used Eqs. (IS4l) in the last step. 

In the general case, using the matching condition in Eq. (IS11|) and the continuum wave functions in Eqs. (I83|) and 
(IS8I) . the inverse scattering matrix can be written as S~^ = B~^A with 


0 

Vt 

\-aVT -iy/T^^ae-'^ 


—ay/T —iy/1 — Tae^'^\ 

VT iy/T^ 

0 

0 -1 / 


(S26) 


/-VT -iy/T^ 

ay/r iay/\ — Te~^^ 
0 1 
V 0 -ae-*^ 


0 \ 
0 1 
ay/T iy/1 — Tae®'^ 
-y/T iy/T^ / 


(S27) 


In the following, we will use the outgoing continuum states. Their wave function is obtained by setting one of the 
outgoing coefficients A™*, (7°“*, , C™* to unity and computing the incoming coefficients via the scattering matrix, 

Eq. (15231) . 


II. TRANSITION RATES 


In this section we calculate the ionization and refill rates induced by a weak ac Zeeman field using second order 
perturbation theory. For this, we first derive the Hamiltonian due to the Zeeman field, Eq. (9) in the main text, in 
the unperturbed basis of the wave functions introduced in Sec. HI Introducing the Bogoliubov transformation, 

(S28) 

A 

inserting it into Eq. (8) in the main text, and symmetrizing the resulting expression, we obtain 

Hz = ^a,A'7a7a' + h.c., (S29) 

AA^ 

where yA,A' is given below Eq. (9) in the main text. Note that yA,A' = —Ea^,a- This allows writing Eq. (jS29p as 

Hz = ^+,-7+7- + E Vv^Efj."tu"tEfj. + 2 E ^Efi.,E'fi'7Efi'lE'fi' | + h-C. (S30) 

Using the definition of ionization and refill rates in the main text, we can calculate them by applying Fermi’s golden 
rule to the Hamiltonian (IS30I) . 

= 27T{fiBBf Y.S{n + E, + E) e{TEE,). (S3I) 

E :f-l 

Here, the upper sign is for an ionization process, whereas the lower sign is for a refill process. The Heaviside function 
appears due to the Fermi-Dirac distributions at zero temperature, ensuring that in a refill process the bound state 
and the continuum state are empty, and in an ionization process the bound state is occupied whereas the continuum 
state is empty. Noting that sign(£’i.) = sign(sin </)), we obtain 

ry^in) = Y l^-;Ti5'sign(sin0).M|' ^ + {\Ey T if')sign(sin0)]. (S32) 

E'>0,fj, 


Using the density of states in the leads p{E)^ we can replace the sum by an integral and obtain Eq. (17) in the main 
text. 

To obtain the matrix elements V^^e^i, we use the general expressions for the bound state wave functions and the 
continuum wave functions, defined above in Eqs. (IHll), dMl) and (IS3L (IS8L respectively. After integration over the 
real space coordinate and reorganization, we obtain 


Vv^E^i = 
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ik 


— ik 


(S33) 
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where 

+ (i?hC”e*^ + D,C^J^e-^^){a + (S34) 

F^ = - -BhAr‘)(ae-*^'' - 1) + - D,C^'^^){1 + ae”*^''). (S35) 

Using the outgoing wave functions as defined above, Eq. (IS33|) can be evaluated numerically to obtain the rates for 
arbitrary transmission. In the following, we consider the two special cases of a transparent and an opaque junction, 
for which we give the analytical expressions. 


A. Transparent junction 


Here, we want to calculate the rates given by Eq. (17) in the main text for T = 1. We have seen, that there are two 
bound states, labelled by R and L, with energies given in Eqs. (IS14I) and (IS15I) . For each of them, we need to evaluate 
(IS33I) . Since the calculation is very similar in both cases, we will only show the explicit calculation for v = Jj. Using 
the coefficients for T = 1, given below Eq. (IS15I1 . and xl = <('/2, we find 


VL,Efj. — 


Kh 


2 ( 1 + a 2 ) 


(a - - C'j)' + l) - l) - 6“*^(l + ae-^i'j 


kl + ik 


kl — ik 


(S36) 


For the transparent junction, there are only two outgoing states, C™* = 1 or = 1. Then, 
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where 


Hi = {a- )(-52*i(ae-*^ + F^) + 53*i(-l + aF^)), 

H2 = -(1 + ae-*^ )(-524(ae-*^ + 6“*^) + 534(-l + 06“*^)). 
Using the unitarity of the scattering matrix, 5i/c5*^. = 5ij, we hnd 
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which after lengthy, but straightforward, algebraic manipulation yields 
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sm • 


(S41) 


where e = E/A. Repeating the calculation for ;/ = R and using the expressions for the bound state energies, Eqs. 
(|S14p and (IS15|1 . we finally obtain 


= 4|e,|(e2-l) 


A(g 2 _e 2 ) 2 (, 2 _, 2 ) 


1 - 


sm 0 es 


(S42) 


where e,^ = E^/A and R = L and L = R. Substituting (IS42|) into (IS32|) and using the energy conservation condition, 
|e| = |U| + je^l for an ionization process and |e| = |U| — \ev\ for a refill process, we find the rates 


= 


(mb-B) 


2 4|ep 


(|U| ± |ei .|)2 - 1 [|U| ± \e„\ +sgn(f 2 ) |ep|sin( 


|U|2(|U|±2|e.|)2[(|U|±|e.|)2-|e,| 


(S43) 
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where we defined 17 = 17/A. 

Near the threshold frequency, |17i{^| = 1 =p |e,y|, the rates grow as \/^, where (517 = |i7| — |17y{^|. At large 
frequencies they decrease as 1/17^. To describe the intermediate frequency regime, we concentrate on the regime 
/>, (517 <C 1 (a similar situation occurs for phases (/ close to tt), where 


pi/R _ 2-\/2\/^ i/R {^bBY \/^(/ 

- - A sn + ^ + “ A 72(^0 +^)2- 


(S44) 


Thus we find that the former reaches its maximum, r__max/[(A^B7?)^/A] = 4/^, at |(517__i„ax| = while the later 

reaches its maximum, r+_niax/[(/rB7?)^/A] = 6-\/3/(/^, at |(517+^max| = (/^/24. 

At larger cj), the maximum is less pronounced and further away from the threshold than for small (p. Note, that for 
(p = 7r/2, since the bound states are degenerate. 


B. Opaque junction 

Here, we evaluate (IS33I) for the opaque junction, T = 0. Using the coefficients derived in Section irm we obtain 

+ (7; 


y-,EtJ. = e 


Av_ 


2(1+ a2) 
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+ ik 


^out _ /^ont picf) 

-(1 + la) + - y (1 - ia) 
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(S45) 


As in the transparent case, there are only two outgoing states, = 1 or = 1. After some algebra and using 
the unitarity of the scattering matrix, we find 




Eti\ 


1 + 


1 + q;^ 1 + k"^ 


(1 — ia^ 


(k_ — ikY 


{S 22 - 524 + e-'^544 - e”‘^542) 


Using the scattering matrix (IS25I1 . we finally obtain 


i2 16(e^ — 1) 

2_^\V-,eA —. 


(S46) 


(S47) 


Since the bound state energy is zero, the rate is identic for refill and ionization processes and reads 


r = 




A 


|17|s 


(S48) 


Near the threshold frequency, |17c| = 1, the rate grows as Vsfi. It reaches its maximum, rmax/[(MBS)^/A] = 
2®/(25-\/5), at |17max| = '/5/2 and decreases as 1/17^ at large frequencies. Note, that the rate does not depend on the 
superconducting phase difference, since the bound state energy is independent of the phase. 


III. STATIONARY OCCUPATIONS 


In this section we give the analytical expressions for the stationary occupations p®/ of the state |1), that are used 
for the calculation of the supercurrent. These stationary occupations are obtained from the master equation, Eq. (II) 
in the main manuscript, by setting p = 0. As the current in the main text, they are most conveniently expressed in 
the form 

Pil =Pu + p2 _p^2(5a;)^ - Pu)- (S49) 

The occupations far from resonance (p^) are given as 
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(S50) 
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whereas the occupations at the resonance (p® ) take the form 

1 2 r / pR 
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The expressions for the current given in Eqs. (12) and (13) in the main manuscript 
stationary occupations, Eqs. (IS51|) - (IS54|) . into Eq. (10) in the main manuscript. 
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are obtained by inserting the 
























